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PERFORMANCE  OF  THE  TRANSITION  ANALYSIS 
PROGRAM  SYSTEM  ON  THE  ADVANCED 
SCIENTIFIC  COMPUTER 

INTRODUCTION 

During  the  past  three  years  the  Naval  Sea  Systems  Command  (NAVSEA 
03)  has  sponsored  the  development  by  McDonnell  Douglas  Corporation  of  a 
large  and  quite  general  computer  code  for  the  analysis  of  flows  around 
submerged  bodies.  It  includes  the  potential  flow  analysis,  boundary 
layer  solution,  and  prediction  of  the  growth  rate  of  disturbances  in 
the  laminar  boundary  layer  (from  which  the  location  of  the  laminar-to- 
turbulent  transition  can  be  predicted).  Users  of  the  code,  called  the 
Transition  Analysis  Program  System  (TAPS),  have  reported  that  it  is 
undesirably  expensive  to  run  on  conventional  computers.  Therefore, 
work  has  been  initiated  at  the  Naval  Research  Laboratory  (NRL)  to 
investigate  possible  cost  reductions  by  using  its  Texas  Instruments 
Advanced  Scientific  Computer  (ASC),  a "pipeline"  machine.  This  work  is 
potentially  beneficial  to  fluid  dynamicists  at  a number  of  Department  of 
Defense  installations  by  virtue  of  their  access  to  the  ASC  through  the 
Defense  Advanced  Research  Projects  Agency  Network  and  anticipated 
access  through  the  Navy  Laboratory  Computer  Network  (NALCON). 

The  present  report  summarizes  work  done  to  date  with  TAPS  on  the 

ASC.  Each  of  the  106  subroutines  of  TAPS  has  been  compiled  with  the 
Note:  Manuscript  submitted  June  26, 1978. 
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vector  (NX)  compiler;  and  the  axisymmetrlc  potential  flow,  laminar 
boundary  layer,  and  transition  portions  of  the  code  have  been  tested 
with  sample  problems.  Comparisons  of  the  results,  run  times,  and  costs 
with  those  obtained  at  other  computational  facilities  have  been  made 
where  possible,  allowing  an  assessment  of  the  advantages  of  using  TAPS 
on  the  ASC. 

REVISIONS  OF  THE  CODE  FOR  ASC  USAGE 

An  IBM-compatible  version  of  TAPS  was  provided  by  NAVSEA.  The 
following  changes  were  made  for  ASC  usage:  (a)  a number  of  constants 
which  were  undefined  in  the  original  code  such  as  NTYPE,  were  either 
eliminated  or  set  equal  to  zero;  (b)  calls  to  IBM  subroutines,  such  as 
ERRSET,  were  eliminated;  (c)  the  write  statements  to  Tape  17  in  sub- 
routine FLOWS  were  changed  from  formatted  to  unformatted  to  be  consis- 
tent with  the  remainder  of  the  program;  (d)  the  READ  and  WRITE  state- 
ments to  Tape  23  in  subroutines  STAB,  OTPT,  and  INP2  were  changed  from 
formatted  to  unformatted. 

Subroutines  were  compiled  at  the  K level  of  optimization  except 
where  difficulties  in  compilation  or  execution  ensued.  Subroutine  OSI 
had  to  be  compiled  at  the  I level,  for  example,  because  higher  level 
compilations  introduced  errors  by  vectorizing  short  DO  loops  with  com- 
plex arithmetic  operations.  It  must  be  emphasized  that  no  changes 
have  been  made  to  date  in  TAPS  to  facilitate  vectorization.  Thus, 
while  a significant  amount  of  vectorization  was  done  by  the  compiler, 
much  more  would  be  possible  by  altering  portions  of  subroutines  to 
facilitate  this  operation. 
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SAMPLE  PROBLEM  ONE 

This  problem  was  used  by  McDonnell  Douglas  to  test  TAPS.  It  con- 
sists of  the  analysis  of  the  flow  around  an  axisymmetric  body  of  approx- 
imately 4.3  m length.  The  maximum  diameter  of  approximately  0.52  m is 
2.7  m from  the  nose,  and  a weak  positive  pressure  gradient  and  a heated 
wall  are  present  over  most  of  the  forebody.  The  program  was  used  to 
generate  the  potential  flow  solution  around  the  body,  predict  the 
laminar  boundary  layer  characteristics  forward  of  laminar  separation, 
and  determine  the  spacial  amplification  ratio  of  selected  disturbances 
to  the  laminar  boundary  layer  at  a freestream  Reynolds  number  per  meter 
of  1.125  x 107. 

The  potential  flow  and  boundary  layer  solutions  of  this  problem 
generated  by  the  ASC  were  essentially  identical  to  those  obtained  pre- 
viously at  other  computational  facilities.  Excellent  agreement  was 
also  realized  in  the  values  of  spacial  amplification  ratios  computed 
on  the  ASC  and  elsewhere.  Table  1 shows  the  amplification  ratios 
computed  for  particular  stations  (axial  locations  on  the  body)  and 
normalized  disturbance  frequencies  on  three  computers.  Similar 
correspondence  in  the  computed  ratios  has  been  verified  at  other 
stations  and  frequencies. 

Detailed  run  time  and  cost  comparisons  have  been  made  between  the 
ASC  and  the  CDC  6600  system  at  the  David  W.  Taylor  Naval  Ship  Research 
and  Development  Center.  Table  2 shows  shorter  run  times  on  the  ASC  and 
Table  3 shows  greatly  reduced  run  costs.  Results  have  been  shown  with 
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both  formatted  and  unformatted  Tape  23  runs  for  the  ASC,  while  CDC  6600 
time  and  cost  Information  are  available  only  with  a formatted  READ  and 
WRITE  to  Tape  23.  The  times  and  costs  for  the  stability  (amplification 
ratio)  analysis  include  five  program- selected  frequencies  and  one  user- 
specified  frequency.  Times  and  costs  refer  only  to  the  central  pro- 
cessor of  the  computational  facility,  as  the  requirements  associated 
with  peripheral  equipment  such  as  printers  depend  on  user-specified 
options  of  TAPS. 

SAMPLE  PROBLEM  TWO 

The  ASC  version  of  TAPS  (unformatted  Tape  23)  was  used  to  study 
the  flow  around  an  axisymmetric  body  with  a strong  positive  pressure 
gradient  on  the  forebody.  Potential  flow,  boundary  layer,  and  spacial 
amplification  calculations  were  conducted  at  a total  of  eleven  volu- 
metric Reynolds  numbers  in  the  .7  x 10^  to  8 x 10**  range.  At  each 
Reynolds  number  the  growth  of  disturbances  of  ten  frequencies  was 
computed  to  establish  the  maximum  amplification  ratio.  The  resulting 
plot  of  maximum  amplification  ratio  as  a function  of  Reynolds  number, 
given  in  Figure  1,  is  important  for  estimating  the  velocity  or  volu- 
metric Reynolds  number  at  which  transition  moves  forward  from  laminar 
separation.  At  this  point  a rapid  increase  in  the  drag  coefficient 
with  Reynolds  number  begins.  Assuming  the  body  surface  to  be  smooth 
this  critical  volumetric  Reynolds  number  in  Figure  1 is  approximately 
4.25  x 10**,  at  which  the  amplification  ratio  is  e9.  The  central 
processor  times  and  costs  for  the  analysis  at  each  volumetric  Reynolds 
number  varied  from  2.0  to  2.2  minutes  and  $29  to  $34,  respectively. 
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(VOLUMETRIC  REYNOLDS  NUMBER  ) x 1CT6 

Figure  1 


DISCUSSION  AND  CONCLUSIONS 


The  work  reported  herein  has  established  the  TAPS  can  be  used 
much  more  economically  on  the  ASC  than  on  other  Navy  computational 
facilitief  Substantial  further  improvement  may  be  realizable  by 
rewriting  selected  subroutines  to  increase  the  amount  of  vectorization. 
The  performance  of  the  two  dimensional  potential  flow  and  turbulent 
boundary  layer  portions  of  TAPS  on  the  ASC  remain  to  be  established. 
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